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Abstract 

TeV-mass dark matter charged under a new GeV-scale gauge force can explain electronic cosmic- 
ray anomalies. We propose that the CoGeNT and DAMA direct detection experiments are ob- 
serving scattering of light stable states — "GeV-Matter" — that are charged under this force and 
constitute a small fraction of the dark matter halo. Dark higgsinos in a supersymmetric dark sector 
are natural candidates for GeV-Matter that scatter off protons with a universal cross-section of 
5 x 10 -38 cm 2 and can naturally be split by 10-30 keV so that their dominant interaction with 
protons is down-scattering. As an example, down-scattering of an 0(5) GeV dark higgsino can 
simultaneously explain the spectra observed by both CoGeNT and DAMA. The event rates in 
these experiments correspond to a GeV-Matter abundance of 0.2 — 1% of the halo mass density. 
This abundance can arise directly from thermal freeze-out at weak coupling, or from the late decay 
of an unstable TeV-scale WIMP. Our proposal can be tested by searches for exotics in the BaBar 
and Belle datasets. 
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I. INTRODUCTION 



For decades the WIMP — a particle with TeV-scale mass that freezes out through weak- 
strength couplings — has been the canonical dark matter candidate. The electronic cosmic 
ray excesses observed by PAMELA p], Fermi [2], and HESS (3JH], along with the discovery 
of the WMAP and Fermi "haze" [5H9] , have both bolstered and complicated the WIMP 
picture. Together these experiments suggest that WIMP dark matter either annihilates or 
decays into lepton-rich final states (see e.g. [10J), with a rate 10-1000 times larger than 
expected from s-wave annihilation. Dark matter charged under a new GeV-scale U(l)u 
'dark force' that kinetically mixes with hypercharge can explain these results through either 
Sommerfeld-enhanced annihilation pTHT5] or decay [H] into light dark-sector states, which 
are in turn forced by kinematics to decay into leptons. Such a GeV-scale sector is naturally 
realized in models with supersymmetry-breaking near the weak scale, where D-term mixing 
dynamically generates a GeV-scale mass for the U(l)r> gauge boson A' [171421] (summarized 
in Section |ll|) . 

The spontaneous breaking of this U(l)r> gauge symmetry requires a higgs sector, and 
this sector can naturally contain stable GeV-scale states such as dark higgsinos. These 
and other stable "GeV-Matter" particles charged under U(1)d interact with protons with a 
cross-section 

a XiP «5x 10~ 38 cm 2 (1) 

that is completely fixed by Standard Model couplings [18]. Because this cross section 
is so large, even a small abundance of GeV-Matter can be observed at direct detection 
experiments. In fact, we should expect GeV-Matter states to comprise only a small fraction 
of the dark matter halo because their annihilation cross-section is large (~ a 2 / ra Gev)' 

Two experiments have reported anomalous events that may arise from scattering of a 
light dark matter species. DAMA J22H2S] has reported a 9a annual modulation signal at 
2-5 keVee in their Nal crystal detector. Early this year, the CoGeNT collaboration reported 
about 100 events from an unknown source in a low-threshold, high-resolution Ge detector 
[26] ; these events are consistent with light dark matter scattering with a mass and cross- 
section similar to those needed to explain the DAMA modulation. As shown in Figure [I] (see 
Section [V] for details), both signals can be explained by 3-12 GeV dark matter scattering 
with a cross-section a c g ~ (1 — 5) x 10 -40 cm 2 , or by a fraction of the halo scattering with 
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a larger cross-section (see also [26H28] ). If the events reported by these experiments arise 
from light dark matter scattering, they provide further evidence for a new sector at the GeV 
scale. 

In light of the cross section ([!]), dark sector higgsinos at 3-12 GeV naturally explain the 
CoGeNT and DAMA excesses if they comprise a fraction 

/n g ht = ^ « (2-10) x 1(T 3 (2) 
Pdm 

of the halo mass density. We propose two natural origins for an abundance of this size: 
Thermal freeze-out generates an abundance (assuming s-wave annihilation) 

f q v in -3 f light W 1/150 V 

/ii g ht~3xl0 j (—J , (3) 

which yields the desired GeV-Matter relic density for small couplings, as we discuss 
in Section IIHI 

Late Decays of TeV-mass WIMP-sector states generate a density 

/hght - m ' ight Phcavy - few x 10- 3 (4) 

^heavy m hea,vy PDM 

of GeV-Matter that can easily dominate over its thermal abundance. This sharp 
prediction is insensitive to the details of the decay process, and agrees with the halo 
fraction required to explain the CoGeNT and DAMA signals! The decaying TeV-mass 
particle can be the scalar superpartner of a stable fermionic WIMP, as discussed in 
Section [TVl 

We therefore propose a new unified framework for interpreting both the direct detection and 
cosmic ray anomalies from the same dynamics. 

£7(1) ^-breaking mass splittings of 0(few keV — few GeV) for both WIMPs and GeV- 
Matter arise from generic interactions with TeV-scale states. In the presence of such mass 
splittings, dark matter scatters through inelastic up- and down-scattering, with elastic pro- 
cesses highly suppressed. These splittings have two important effects: when the GeV-Matter 
states are split by O(10 keV), direct detection is dominated by down-scattering, which im- 
proves the agreement between DAMA and CoGeNT and permits lower masses of the light 
dark matter than in the elastic discussed in Section M and more generally in [29J . 
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FIG. 1: Left: Regions of scattering cross-section per nucleon times GeV-Matter mass fraction 
in the halo versus mass favored by CoGeNT (filled regions), DAMA with channeling (solid- line 
contours), and DAMA without channeling (dashed contours) for light dark matter scattering co- 
herently off protons. We show favored regions (from right to left) for elastic scattering (green) and 
inelastic splittings |<5| = 15 keV (orange), 25 keV (red), and 35 keV (blue), with equal fractions of 
up- and down-scattering. The contours correspond to "1.6<r," neglecting systematic uncertainties. 
The black dot and upside-down triangle correspond to benchmark points we consider in Figure [2} 
Right: A close-up of the |<5| = 25 keV allowed region with "3<r" contours added and constraints 
from XENON10 (black lines) and CDMS-Si (gray lines). Thick (thin) lines denote 5.0 (2.3) events 
expected, again neglecting systematics. More details are given in E jVj 

A larger > MeV splitting of the TeV-scale WIMP allows it to explain the INTEGRAL 511 
keV excess and renders it practically invisible to direct detection experiments [TT| 150]. 

The success of these models is particularly striking in light of the difficulty of generating 
such high direct detection rates through Standard Model processes, even when a large abun- 
dance of light dark matter is generated through an asymmetry [31] . Dark matter coupled to 
the Z boson can explain the reported rate, but is quite constrained by LEP measurements 
of the Z invisible width [32] ■ Standard Model Higgs exchange, in contrast, leads to typical 
cross-sections < 10~ 43 cm 2 (see e.g. [33]) because the Higgs couples only weakly to nuclei. 
Explanations of the CoGeNT signal through scalar interactions posit a new light scalar that 
mixes significantly with the Higgs, but remains relatively light through a percent-level fine- 
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tuning [27J|28]. I n contrast, the low mass of our dark matter and mediator arise from simple 
dynamics. 

Implications for S-factories 

Analysis of existing low energy collider data is crucially important for testing the GeV- 
Matter scenarios in this paper (see e.g. p£fH55] ). The strongest constraint on U(1)d gauge 
boson decays to Standard Model matter is from a BaBar search for a resonance in 7/i + /i~ 
in T(3s) data jlQ], which is sensitive to the "radiative return" process e + e _ — > jA'. The 
parameter space consistent with the GeV-Matter explanation of the CoGeNT and DAMA 
excesses is just below the limit from this search, but can be tested by a search in the 
full Belle and BaBar datasets, which together comprise about a factor of 30 higher 
statistics than the T(3s) data. The "higgs'-strahlung" process |35J may also be accessible 
at B- factories, leading to a striking six-lepton final state. Searches at J/ip- and ^-factories 
can also explore the low-mass parameter space [36j E7J HU 02] . 

Down-Scattering and Lighter Dark Matter 

Three flaws have been noted with elastically scattering light dark matter as an expla- 
nation of the DAMA annual modulation signal: the predicted modulation spectrum grows 
exponentially at low energies while the DAMA spectrum appears to fall near threshold [43J . 
Moreover, the mass range favored by both the DAMA and CoGeNT spectra, ~ 10 GeV, is 
disfavored by XENON10 HMS] and especially CDMS Silicon H7J data (see [27]). Finally, 
the scattering rates required to explain the CoGeNT and DAMA signals differ by about a 
factor of 5 (though this final discrepancy depends on the precise channeling properties of 
Iodine [STJ HHH50]). 

However, there is no reason to expect the light dark matter scattering to be purely elastic. 
Mass splittings of order 10-30 keV for GeV-Matter are readily generated by TeV-suppressed 
higher- dimension operators. 

Because the dark sector kinetically decouples from the Standard Model at temperatures 
high compared to this splitting, both states will be equally populated and the heavy state 
will be cosmologically long-lived [5H453] . The excited states interact with nuclei through 
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down-scattering, which pushes the nuclear recoil spectrum to higher energies for a given 
dark matter mass. As we discuss in £|V] this effect improves the agreement of light dark 
matter with DAMA modulation data (relative to the elastic case), brings the CoGeNT 
and DAMA expected cross-sections together, and is consistent with the XENON10 null 
search. The residual tension of this scenario with the 5-tower CDMS Silicon analysis |47] is 
extremely sensitive to the experimental selection efficiency near threshold and small energy- 
scale systematics. There is also tension with the DAMA unmodulated spectrum in [24]. 
A general analysis of down-scattering light dark matter [29] also finds a lower-mass region 
consistent with DAMA and CDMS Silicon results. 



II. AN APPROXIMATELY SUPERSYMMETRIC DARK SECTOR 

Our starting point is a simple extension ( [T8Tl2~T] ) of the supersymmetric SU(3) x SU (2)^x 
U(1)y Standard Model that dynamically generates a mass scale of O(GeV) for the "dark 
force." The dark sector consists of two chiral superfields h u and h d that are oppositely 
charged (±1) under a new U(1)d gauge group, but neutral under the Standard Model gauge 
groups. No Standard Model fields are directly charged under U{l)u. We also add a gauge 
singlet S, and consider the following superpotential for the GeV-scale sector: 

W = \Sh u h d - \w 2 D - € -W«W aD . (5) 

Here, Wd and Wy are the U{1)d and hypercharge vector superfields, and the last term is 
the supersymmetric version of the usual gauge kinetic mixing, 

£ = \ Fy^Fg. (6) 

Besides inducing the gauge kinetic mixing @ and gaugino kinetic mixing, the last term also 
induces a fixed e-suppressed potential coupling of Standard Model and dark-sector scalar 
fields through the mixing of Standard Model and dark-sector Z)-terms. After electroweak 
symmetry breaking, this mixing generates a dark-sector scalar potential 

V = A 2 (\sh u \ 2 + \Sh d \ 2 + \h u h d \ 2 ) + y(N 2 - M 2 + C 2 ) 2 , (7) 



that includes a negative mass-squared for the scalar h d (or h u if ( 2 < 0), where 

egy | cos 2(3 \ 



C = ^rT^ l (8) 
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is the hypercharge D-term vev. 

The stable vacuum is one where U(1)d is broken but supersymmetry is approximately 
preserved (supersymmetry-breaking corrections, discussed below, are higher-order in e). The 
light fields are grouped into massive gauge and chiral supermultiplets, with masses 

rriA' = m \ = V2~9d( and m h = m~ h = X(. (9) 

The uneaten real part of h d and h d (which gives a Dirac mass to the gaugino A) are in the 
massive gauge supermultiplet, while the h u and S scalars and their superpartners h u and S 
marry into the massive chiral supermultiplet. The A' mediates scattering of particles with 
£/(l)£>-charge Qd off Standard Model matter with a non-relativistic scattering cross- 
section almost completely determined by Standard Model parameters, 

_ ^ p a D ac 2 w ^ _ (Wnqja^)^ f 2q d \ 2 

a ™ ~ 1&KQd < " (cos 2/3)^ ~ 5 X 10 {^2(3 ) Cm ' 

with all dependence on e and the dark gauge coupling go canceling and only mild dependence 
on the ratio tan/3 of Standard Model Higgs vevs and on the dark matter mass (we chose 
m x = 5 GeV in the numerical estimate above). 

In addition to these light fields, we assume throughout that the dominant WIMP compo- 
nent of dark matter arises from a vector-like pair of chiral superfields $ Ui( j with charge ±q, 
so that the full renormalizable superpotential is 

Wtot = *Sh u h d - \w 2 D - € -W?W aD + M$A, (11) 

with M at the TeV-scale (for example, M and the /x-term of the MSSM Higgs fields may 
have the same origin). The stability of $ and absence of other couplings to the light fields 
can easily be guaranteed by global symmetries, in particular including a Z 2 under which the 
$> U: d fields are odd. Alternately, the stability of $ is automatic if its charge q is non-integer. 



Supersymmetry-Breaking Effects 

Supersymmetry-breaking effects are crucial for dark sector physics; we will always take 
fermionic states as our dark matter candidates, so we will insist that they are stable. In the 
light GeV-Matter sector, some scalar density is innocuous but it is important that the $ 
scalars decay. 
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We begin by discussing the supersymmetry-breaking effects that arise inevitably from 
gauge mediation through the kinetic mixing and from anomaly mediation, and then we will 
consider further effects that can arise through gravity or direct mediation to the dark sector. 

In models of gauge mediation, all scalars charged under U(1)d receive soft masses medi- 
ated from the messengers through kinetic mixing of order [201 G3] 

~ 2 2 a D 2 /-i o\ 

m GMSB ~ e — m Ec , (12) 
ay 

where is the mass of the right-handed selectron in the Minimal Supersymmetric Stan- 
dard Model (MSSM). A negative supersymmetry breaking mass-squared for S is generated 
by RG evolution [20], although it is suppressed by an additional loop factor times a po- 
tentially large logarithm. The i?-breaking Ah u hdS and B^ u ^d terms are suppressed by 
0(ar>e 2 ) relative to the MSSM bino mass, wig. A Majorana soft mass for the dark gaugino 
is given by e 2 m§, while the MSSM bino and dark gaugino have a mixing mass of em^. In the 
absence of U(1)d symmetry-breaking, these would preserve a zero mass eigenvalue. When 
U(1)d is broken they only induce a tiny splitting e 2 m|/m^ of the Dirac gauginos into two 
Majorana states (approximately (A ± h</)/y/2), one heavier and the other lighter than the 
A'. Both the lighter and heavier of these states can decay through kinetic mixing with the 
MSSM bino to 7G, with a lifetime [20] 

1 



in , / y/F \ /lGeV\ 5 /l0- 3 \ 2 

10 s to^v UrJ {—) (13) 



as long as m~ x > m 3 / 2 - Constraints from BBN |55j and from EGRET searches for 7-ray lines 
from the galactic center require [56] either very low y/F < 10 5 GeV for prompt gaugino 
decays or rather high y/F ~ 10 10 GeV, so that the gravitino is heavier than the gaugino 
and the decay is forbidden. In the latter case, the dark sector soft masses inevitably receive 
anomaly-mediated contributions [53, |58] of order 

OLD 

?"AMSB ~ ^—"^3/2 ~ 1 — 100 MeV, (14) 

which are comparable to the gauge-mediated contributions to the scalar soft masses and 
dominate the gaugino soft mass. 

If F is low, dark matter annihilation or decay to A's at a rate compatible with the 
FERMI e + e~ excess could potentially generate a gamma-ray signal visible with existing 
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or upcoming satellite data. This is because supersymmetric couplings will also allow dark 
matter annihilation/decay to A, whose subsequent decay produces a photon. This signal 
is fairly model-dependent - - it is reduced if dark matter annihilates or decays through 
supersymmetry breaking couplings, or if the gaugino decays primarily through i?-parity- 
violating interactions. Nonetheless, it would be interesting to investigate these signals in 
future work. 

Although the effects we considered above are irreducible, it is plausible to include addi- 
tional supersymmetry breaking operators such as 



where X is a supersymmetry-breaking-spurion with an F x component vev. We can think of 
these terms as originating from 'gravity mediation' with M* ~ M ph but it is equally natural 
for M* to be much smaller. These operators can contribute to the Majorana gaugino mass 
and also raise the scalar higgs sector (S and h u ) soft masses. In particular, the scalar S 
can obtain a net positive supersymmetry-breaking-mass, forbidding the decay of the lightest 
higgsino to S and a gravitino. 

Stable States and Their Properties 

We will always take our GeV-Matter (the light stable states) to be fermionic, so in this 
section we will be interested in the properties of these fermions and in the mechanisms 
through which we avoid stable scalars. 

The Dirac gaugino can be stabilized on cosmological timescales if the supersymmetry- 
breaking F is large, so that either the lifetime (oc F 2 ) of the gaugino is very large or 
m 3 / 2 > m~ x (which requires a/F > 0.5 x 10 10 GeV) so that it is kinemtically forbidden to 
decay into the gravitino. Nonetheless, it is not a good candidate to explain the CoGeNT and 
DAMA signals in our scenario. For such high F, the gaugino anomaly-mediated soft mass 
~ 10 — 100 MeV is sufficiently large that the heavier gaugino decays rapidly to e + e~ and 
the lighter gaugino. Though the lighter component can comprise a non-negligible fraction 
of the dark matter, it can only scatter elastically (because the excited state is depopulated 
and up-scattering is kinematically forbidden) with a cross-section that is velocity-suppressed 
and thus too small to explain the direct detection signals. 




(15) 
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The higgs supermultiplet thus contains the only viable GeV-Matter candidates within the 
model ([5]) that can explain CoGeNT and DAMA. We assume that the higgs supermultiplet 
is stabilized on cosmological timescales by a discrete or continuous global symmetry. We 
are primarily interested in higgsino GeV-Matter {{h u ± «9)/v2), which readily acquires a 
small inelastic splitting as we discuss below. The higgsinos are stable under either of two 
conditions: either S must receive a positive soft-mass contribution that makes it heavier 
than the higgsino, or F must be sufficiently large that decays h — > S + G are kinematically 
forbidden. 

If F is low and the supersymmetry breaking contribution to the S scalar mass is 
positive (which requires some contribution besides gauge mediation), all scalars can de- 
cay to their fermionic partners and a gravitino with a lifetime given parametrically by 



t ~ 167rF 2 /(m^Am), and likewise gauginos can decay to 7G according to (13). The hig- 
gsino is then the only stable relic in the light dark sector. In this case, the scalar $ also 
typically decays to $ + G. 

Large F also leads to a stable higgsino if the gravitino mass is large enough to kinemat- 
ically forbid higgsino decays to S + G or S + ha- In this scenario, one or both scalars in 
the higgs multiplet may be stable, but they are less phenomenologically important than the 
higgsinos. In particular, S is neutral under U(1)d and therefore only has a very small direct 
detection cross-section induced by its mass mixing with h u . Even at high F, a heavier h u 
scalar can typically decay through the same mass mixing to S and an off-shell A'. Even if the 
h u is cosmologically stable, it typically has no inelastic splitting, making it less observable in 
direct detection than the higgsinos which, as we will discuss below, typically down-scatter. 

A final possibility that we will not explore further in this work is that the higgsinos are 
actually lighter than the gaugino and gravitino states, so that the higgsinos are stabilized 
by R-parity. 

Although supersymmetry-breaking effects leave the higgsino and $ as exactly degenerate 
Dirac fermions (in contrast with the gaugino), interactions with 1-10 TeV states can lead to 
splittings by O(10 keV — 1 MeV) when these heavy states are integrated out, for example, 
through the higher-dimension superpotential operator 

< 16 > 

Elastic scattering processes will be velocity suppressed and therefore negligible in all cases, 
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so the higgsinos will only be visible to direct detection experiments through inelastic down- 
scattering. 

In many models (see Appendix [B]) , the same interactions generate an 0( MeV) or larger 
splitting for the heavy dark matter fermions <& u ,d and a smaller (9(10 keV) splitting for the 
higgsinos; the simplest example is to take q = 1/2 and include the marginal superpotential 
operator h u $ 2 d . For splittings SM > 1 MeV, the excited state of the TeV-scale dark matter 
decays to the lower-mass state and an e + e~ pair, depopulating the initial state of the down- 
scattering process, and up-scattering off nuclei is kinematically forbidden. In this case, 
the leading direct detection cross-section would be elastic, but it is suppressed by both 
(^) ~ 1(T 12 - 1(T 6 and by a factor of the non-relativistic velocity squared ~ 10 6 , making 
it completely negligible. For splittings ~ 1 MeV the <& u ,d self-interactions may explain the 
INTEGRAL 511 keV excess [EE]. 

We note that it is less generic in these sectors to generate inelastic splittings for the scalar 
$ than for the $ fermion, so we insist that the scalar $ decay. At low F the gauge-mediated 
splittings are sufficient to allow decays $—)■$ + G. For the higher- F scenario and heavier 
gravitinos, it is plausible that the $ couples more directly to SUSY-breaking than do the 
fields in the low-mass dark sector (for example, if its /i-term is related to SUSY-breaking), 
so that decays to the gravitino are still kinematically allowed. Alternately, in Section [TV] we 
consider interactions with an intermediate-scale supersymmetry-breaking sector that lead 
to $ scalar decays but keep the fermion $ cosmologically stable. 

III. GEV-MATTER FROM THERMAL FREEZE-OUT 

Given the proximity of the CoGeNT and DAMA mass scales to the GeV-scale, a natural 
possibility is that the higgsinos of the dark sector are themselves stable and give rise to 
observable direct detection signals. We focus on A > V2gD, so that the higgs supermulti- 
plet is heavier than the gauge supermultiplet, and we will assume that global symmetries 
guarantee that this multiplet is stable on time-scales of order the age of the universe. Note 
that the gaugino states could be cosmologically stable, but they are not a viable candidate 
to explain the CoGeNT and DAMA signals as we discussed in Section [TTJ The higgsino relic 
abundance agrees with ^ in a small-coupling limit. For larger couplings the thermal hig- 
gsino abundance is small but, as we will see in the next section, the late decay of a TeV-mass 
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WIMP gives rise to a relic density of higgsinos of the appropriate size. 

We will also discuss how an additional light vector-like species (not part of the dark gauge 
or higgs sector) is a viable GeV-matter candidate if it is charged under U(1)d- 



Dark Higgsino Freeze-Out 

Higgsino annihilation during freeze-out is dominated by annihilation to the real scalar hd 
and a longitudinal A', with cross-section 

a h = a 2~~ I 17 ) 

h 

where a\ = X 2 /(Att). This dominates over all other annihilation modes, which are either 
p-wave suppressed or suppressed by g 2 /\ 2 , leading to an abundance 

- 3 (^) 2 (aiv) 2 - W> 

where ot\ = X 2 /(Att), and we have chosen the baseline m~ h = 6 GeV near the CoGeNT 



and DAMA best-fit regions for 25keV inelastic splitting (see Section [V]). Using (10) with 
qd = 1/2, this gives rise to an effective direct-detection cross-section 

a cS = f- h a- hp ~ 1.5 x 10- 40 cm 2 (19) 

x 



/ l V/l/lSOV/ y 
Vcos2/3y V "a J V6GeV/ 



Normalizing this to the a e g derived from DAMA and CoGeNT data determines the dark- 
sector higgs vev within this setup (at around 20 GeV), and the gauge coupling as a function 
of the A' mass: 

aD = 2 a "{^) " 10 (iGiv) (eGev) ( l.5xl0-«c m ^ (20) 

near the CoGeNT/DAMA best fit. A separation of scales between the A' and higgsino 
masses, as suggested by the lepton-rich cosmic ray signals, implies a rather small dark gauge 
coupling in this scenario. This in turn can be translated into a prediction for e using ^ 
and g, 

5 * (ilv) W^r^^f • (21) 
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If the A' decays directly to e + e , models that produce <r e g as a thermal relic are very near 



the sensitivity limit of existing 5-factory searches discussed in Section VI 

For larger gauge couplings, the relic abundance of the higgsinos from thermal freeze-out 
is too small to explain the CoGeNT and DAMA signals. However, as we will discuss in 



Section IV, a relic abundance of the desired size can easily be generated by other means. 



Vector-Like GeV-Matter 

Now let us consider a slightly less minimal model that exhibits the parametrics of equation 
(|3j) without a direct link between the GeV-Matter mass and the A' mass. Consider a light 
vector-like species x± with superpotential W D [i x x+X- an d H x ~ 5-10 GeV. Like <& u ,d, 
the x± states may be stabilized by a discrete symmetry such as a Z2 where x± ~~ >* ~X±- A 
higher-dimension coupling of x to the higgs-sector field hd such as 

\x\hl (22) 

will lead to a ~ 10 keV splitting for the fermionic light dark matter states x± when (hd) ~ 
0.1 — 1 GeV and A > 1 — 10 TeV. Note that the scalar states in the x mult ip lets can decay 
to x + C or x + A pairs depending on the gravitino and gauge supermultiplet masses. 

In a minimal scenario, the x re li c abundance is set by the annihilation cross-section to 
lighter higgs and gauge- mult iplet states, which we parametrize as 

a* ~ ^ 23 

x 2m 2 x v V ; 

This gives rise to a relic density 

h = ^^ 6xW -*(iimY(j^)\ (24) 

x Om V <*d J V6GeV/ v ' 



The x scatter in direct detection experiments with total cross-section (10) with qd = 1 
giving rise to an effective cross-section 

^dm 

/ 1 \ 2 /1 /inn\ 2 / 



^cfr = fx7^- a k P - 1-2 x 10- 40 cm 2 (25) 

2 2 



x 



1 V /1/lOOV ( m- n 1 



cos 2/3 J \ a.£) ) \6GeV/ 
The same parametrics applies to the heavy (again, assuming that annihilation through 
gauge interactions dominates), albeit with a higher in the freeze-out calculation. For 
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example, an = 1/100 used above would lead to Q$ = ft dm for /i$ = 550 GeV. The above 
model is similar to that discussed in [59J, but we consider a different hierarchy of masses 
and couplings. 



IV. GEV-MATTER FROM THE DECAY OF HEAVY WIMPS 

In this section, we will present a simple mechanism that can simultaneously explain 
the cosmic ray data and the direct detection signals. Our mechanism populates the GeV- 
Matter states via the late decay of a TeV-scale WIMP-sector particle, resulting in a relic 
abundance that precisely fits the CoGeNT and DAMA direct detection rates (with the 



scattering cross-section in Eq. (10)); it explains the cosmic-ray data through the decay of a 
TeV-scale WIMP with a lifetime on the order of 10 26 s. Both of these effects originate from 
the same supersymmetric operator. 

We consider a setup where the fermionic WIMP and its scalar superpartner freeze out 
with comparable relic densities. While the WIMP itself will have a lifetime greater than 
the age of the Universe, its scalar partner will decay more rapidly into the dark higgsino 
state, providing the dominant contribution to the GeV-Matter density, as we now show. 
If we denote the mass of the light and heavy states by and m#, respectively, then 
the relic density of the light and heavy states from freeze-out are pL(f. o .)0c o"L(ann.) oc m 2 L 
and p#(f.o.)oc <JH{xan.) ocmjj, respectively, so that Pl{Lo.)/ Ph{Lo.) ~ m 2 L /m 2 H for roughly equal 
couplings. After one of the heavy states decays into the light state, the light state will inherit 
the number density of the heavy state, Ul ~ nn = PH{f.o.)/ m H oc Tnjj, so that the relative 
relic density from decay Pl/ph ~ ?nL/ m H 3> m 2 L /m 2 H dominates over pi,(f.o.)- Moreover, for 
mi, ~5 GeV and m# ~ 1—2.5 TeV, the stable heavy state can dominate the dark matter relic 
density so Pl/ Pdm~ m L/ m H ~ (2 — 10) x 10~ 3 , which is in precise accord with the DAMA 
and CoGeNT rates (see Eq. p])) if the light state has a direct detection cross-section given 



by (10) 



For a concrete example of this scenario, we again consider the low-energy superpotential 

W = -\w 2 D - e -W^W aD + XSh u h d + M^ d + ^h 2 d h 2 u + 2^2*2 (26) 

where the S, h u , and hd are GeV-scale states and the fermions <& Uy d are the TeV-scale WIMPs. 
The two higher-dimensional operators induce the necessary splittings among the WIMPs and 
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GeV-Matter, so that in direct detection experiments the GeV-Matter can down-scatter and 
the heavy dark matter component is forbidden from scattering. Now let us include the 
Kahler potential terms 

/ dH W {^ V huX^Y + $t e v huX Yi + h.c.) , (27) 

where X and Y are heavy fields in the supersymmetry breaking sector, and V is the U(1)d 
vector field. When X and Y acquire non-zero Fx and Fy -terms, we break two Z 2 symmetries 
under which both superfields $ n and $^ and either X or Y are charged, allowing the 
scalar to decay through mass mixing with the h u state. Any non-zero 5^-term, generated, 
for example, from gauge mediation, will maximally mix the $ u and $^ states, and also allow 
the $rf scalar to decay. Two residual (non-supersymmetric) accidental Z2's under which the 
fermions $ n and $^ and either the X or the Y scalars are charged, are only broken by 
the scalar vevs (X) and (Y), and induces <& Ut d fermion decay through kinetic mixing with 
h u . The fermions are parametrically longer-lived than the scalars, because the higher mass 
dimension of Fx, Fy must be compensated for by the TeV-scale $ mass. 

The dominant <& scalar decay channels are then $ — y h u + A', $ — > h u + A, and especially 
$ — > hd + S due to the fact that A > Qd] we estimate the lifetime as 

- ~ U) ^?vF (28) 

/1/Q0\ / A/T \4 /n.. m7n„ir\8 / a r \ 3 

10 4 S 



A/30A / M* y / 3 x 10 7 GeV y / M V 
V « A ) Vl0 19 GevJ V VF J \2TeV ) 



where we have set Fx = Fy = F. The dominant $ fermion decays are $ — > h u + A', 
$ -)■ /i u + A, $ ->• /it + S*, and $ -)■ h* d + S 1 " 1 " with lifetime 



^1/30^ M 4 

1/30 \ / M* A 4 /(2 x 10 6 GeV) 2 y /2TeV 



T * ~ (—J Hfi (29) 
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10 " ' a A ) Vl0 19 GeV7 V W(^) 7 V M 
It is important to this argument that less suppressed operators involving X and/or Y and 
Kahler terms involving chiral products such as & u hd are forbidden (as in the UV completion 
in Appendix A ) or that one of the scalar vevs (X) , (Y) vanishes. In either of these situations, 



the interactions in (27) are the leading decays. In general, all that is required to populate 
the higgsino states is for scalar $ decay to occur, so fermion $ decays can be rendered 
negligible if (X) < y/F or (Y) < y/F. 
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The scalar $ decays occur well after the dark higgs and gauge sectors freeze out, so that 
the abundance of the GeV-Matter and potentially long-lived (depending on the gravitino 
mass) A, hd are determined by the density of $. Though initially energetic, the $ decay prod- 
ucts redshift and become non-relativistic before matter-radiation equality. Consequently, the 
velocity profile of these decay products in galaxy halos is controlled by gravitational infall 
like that of the WIMP dark matter. A lifetime in the range 10 s < r$ < 10 5 s conservatively 
avoids BBN constraints for Q$h 2 ~ 0.1 and assuming that $ decays never produce baryons; 
if the $ branching fraction to Standard Model leptons is small (0(1O~ 2 ) is readily achieved 
from ratios of couplings) then shorter and longer lifetimes are allowed [55] . 

The fermion decays $ — >■ h u + A', $ — > h u + A and $ — > h* d + produce leptons through 
the decay A -> £+£~, h d ->■ £+£~, and h u -» S + A -> S + £ + t or h u -> S + £ + £~ 
(the latter decay is through an off-shell A if mh u —m s < m^). Consequently, we expect 
that this model can explain the PAMELA and FERMI signals as arising from $ decay for 
(X) ~ (Y) ~ y/Fx ~ y/Fy ~ 10 3 — 10 4 TeV, with a possible higher-energy component of the 
cosmic ray spectrum arising from Sommerfeld-enhanced annihilations (which by itself may 
not suffice to explain the signal [60]). Note that the $ decays also produce high-energy 
gauginos, which can decay into high energy gravitinos and photons, A — )■ 7 + G. If this 
decay is fast, then it may generate a photon signals close to current limits (see e.g. [6TH63] ). 
This constraint is trivially avoided if the gravitino is heavy and/or A is long-lived. (Note 
that if the gaugino is unstable, its lifetime should also be < 10 5 s to avoid constraints from 
photodissociation of nuclei [55, 64J) That dark matter decay induced by higher-dimensional 
operators can produce sizable cosmic ray signals was first pointed out in [65J, and more 
recently considered in [66], [6?] . 

We present a UV completion of the above model (26) and (27) in Appendix [Al 



V. DIRECT DETECTION PHENOMENOLOGY 



In this section, we discuss direct detection signals from light dark matter species at 
DAMA and CoGeNT, and the constraints from XENON-10 and CDMS-Si. Our focus is on 
extending previous analyses |27] to models of light dark matter with an inelastic splitting 
5 < 35 keV and where down-scattering naturally dominates, which (i) changes the recoil and 
modulation spectra, especially at low recoil energies, (ii) allows lower dark matter masses 
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(down to 3-4 GeV), (iii) permits excellent simultaneous fits to DAMA and CoGeNT spectra 
and rates, and (iv) greatly reduces the strength of limits from XENON10, and magnifies the 
sensitivity of CDMS-Si limits to systematic uncertainties. Down-scattering at lower dark 
matter masses has also been proposed as an explanation of the DAMA signal (29] [68] , and 
[29] analyzes the phenomenology of a variety of down-scattering signals in detail. 

A. Kinematics of Inelastic Down-Scattering 

Vector-current interactions of fermions with a small mass splitting are always inelastic. 
When one component of a Dirac fermion ^ = (ipL ^r) of mass m receives a small Majorana 
mass 5, it splits ^ into two components with Majorana masses ~ m ± 5/2, which 
are approximately ipi = (ipL + ^Pr)/V^ an d tp 2 = i^L — 4>r)/V2- Upon rotation of the 
Dirac fermion's gauge interactions to the mass eigenstate basis, all couplings g^A^^id yifij 
are generated. If the initial state's coupling is vector-like, the coefficients ga are suppressed 
by 5/m (as expected because these couplings break the U{1)d gauge symmetry, and 5 
parametrizes the U(1)d breaking felt by the fermion), but if ipL and ipn have different chiral 
couplings, the ga are generically non-zero. An additional suppression of elastic scattering 
exists irrespective of the size of ga, since the A' couplings to matter are very nearly pure 
vector-like (with axial couplings down by (m^/mz) 2 and proportional to the nuclear spin 
rather than its total charge), and since the matrix element for a Majorana fermion scattering 
off a vector current vanishes at zero velocity. Thus, even for the chiral case, elastic scattering 
cross-sections are always suppressed by v 2 . 

Every halo particle in the high-mass state can down-scatter (unlike up-scattering, there 
is no threshold kinetic energy). However, the recoil kinematics of down-scattering differs 
markedly from that of elastic scattering. In the following, we briefly describe the kinematics 
of down-scattering events for a dark matter particle of mass M x scattering off a nucleus 
of mass rriN- If the dark matter is incident at velocity v, the center-of-mass frame has a 
velocity v cm = fiv/rriN relative to the lab frame, where \i is the dark matter-nucleus reduced 
mass. The squared velocity of the recoiling nucleus in the center-of-mass frame, assuming 
the dark matter mass increases by 5 in scattering (i.e. 5 < for down-scattering) is simply 




(30) 
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Therefore, the average lab-frame recoil energy is 

E(v) = ^m N (vf ab ) = ^m N (v 2 cm + u^ tjCm ) = ^ (v 2 + , (31) 
while the range of recoil energies is given by 

^hFf^'-^hF!)' (32) 

Thus, when |5| > /it> 2 , the effect of down-scattering on the nuclear recoil spectrum is dra- 
matic, with mean nuclear recoil energy 

(E R ) « ^ (33) 



and width AE 1 ^ ~ (jj,/m N )^85j2v 2 . 

Because the recoil spectrum is not peaked at zero, experiments with low energy thresh- 
olds relative to \S\fj,/m N are sensitive to a large fraction of scattering events, even if the 
scattering WIMP is quite light (an elastically scattering state of roughly the same mass, 
such as a dark higgs, can have a comparable halo abundance but go un-noticed because 
its detection efficiency is typically much lower). On the other hand, experiments with high 
energy thresholds, which are only sensitive to the tail of the recoil energy distributions, are 
less sensitive to down-scattering than to elastic scattering with the same characteristic scale 
for E R . This is because, in the elastic case, the typical E R and its exponential tail are both 
determined by (fi 2 /m N )v 2 . Down-scattering models achieve the same characteristic E R scale 
with lower dark matter masses, and therefore a sharper exponential fall-off in rate at high 
recoil energies. Therefore, the importance of differences in energy thresholds between exper- 
iments is amplified. Moreover, the bounds from these experiments are extremely sensitive 
to uncertainties in the threshold energy and the detection efficiency near threshold. 

An approximate figure of merit for comparison of thresholds across different target nuclei 
(when rriDM <C ™>n) is the threshold recoil momentum 

= y/2E Rtthresh m N . (34) 

Reduced mass factors further decrease the sensitivity of detectors using the lightest nuclei, 
but for m ~ 5 GeV the effect is small, even for Silicon. The nominal CoGeNT, CDMS- 
Si, and XENON10 thresholds correspond to = 16.0, 19.2, and 32.2 MeV respectively. 
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Thus, the increased threshold-dependence of downscattering dark matter will significantly 
weaken the XENON10 constraint on downscattering models. The impact on nominal limits 
from CMDS-Si is quite mild, as its recoil momentum threshold is within 30% of CoGeNT's. 
However, the CDMS-Si efficiency near threshold (17] is small and varies rapidly. The above 
comparison makes clear that any uncertainty in near-threshold acceptance affects CDMS-Si 
sensitivity dramatically. 



B. Direct Detection Scattering Rate 

The differential rate for spin-independent scattering of a dark matter species x off a target 

containing nuclei of mass is given by 

dR m N f x p dm (f p Z + f n (A - Z)f 

~ T oT/2~Z. « a x,P 



dE R 2/4 m x 

ft 



i 2 

''■'CSC 



"X 



xF 2 (E R ) / d s v J -^ ^, (35) 

where E R is the energy of the recoiling nucleus, \x n is the dark matter-nucleon (proton or 
neutron) reduced mass, pam — 0.3 GeV/ cm 3 is the local dark matter density, m x and f x are 
the mass and fractional abundance p x /pdm of Xi % an d A are the atomic number and atomic 
mass of the target nucleus, f p and f n are the effective coupling strengths of the dark matter 
particle to protons and neutrons, respectively (for kinetic mixing, f n = so the scattering 
rate is proportional to Z 2 ), and 

°» s udh* f da -\ M r (36) 

is the effective scattering cross-section of the dark matter particle off protons (this equals 
the physical scattering cross-section in the elastic limit, but differs for inelastic scattering 
by a ratio of outgoing to incoming velocities in the CM frame). The parametrization of the 
form factor F(Er), dark matter velocity distribution f(v — v r ), and ^-dependent velocity 



threshold v m i n in the last two terms of (35) are described in Appendix [Cj 

Throughout this section, we treat <J X)P f x as the parameter to be fit or constrained by 
event rates in direct detection experiments. For the inelastic scenarios, f x refers to the total 
X mass fraction, which we assume to be equally divided between the excited and ground 
states, so that the fraction of the halo that can down-scatter is / x /2, while the fraction that 
can up-scatter is f x /2 times a Boltzmann factor, and we add the two contributions. 
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C. The Data 



CoGeNT CoGeNT is an ultra-low-noise germanium detector with a fiducial mass of 
330 g that operated for 56 days in the Soudan Underground Laboratory with an ionization 
energy threshold of 0.4 keVee (keV electron equivalent, i.e. ionization energy). Between 
0.4 keVee and 3.2 keVee, the detector recorded an approximately uniform background rate, 
sharp cosmogenic peaks, and about 100 events of unknown origin between 0.4 and 1.0 
keVee, peaked near threshold [26J. Background explanations of the latter component are 
not excluded, but its spectrum is sufficiently similar to that expected from light dark matter 
scattering off germanium to warrant further study [2"6T[2"8| I3"2"| 1ST?] . 

The ionization energy E ee (in keVee) is related to the nuclear recoil of energy Er by a 
quenching factor. Following [70J, we use a fit to the Lindhard model with k = 0.2, 

M' 12 - (37) 

Assuming this quenching factor, the threshold energy is 0.4 keV ee = 1.9 keV r and = 16.0 
MeV. 

The proximity of the putative signal to the detector's noise threshold and the similarity of 
its shape to that of a hypothetical exponential background pose a challenge for interpretation 
of the CoGeNT dataset. We do not fit to the first bin (0.4-0.45 keVee) on the grounds 
that it is within la energy resolution from the noise threshold and may be significantly 
contaminated. We fit three contributions to the remaining bins between 0.45 and 3.2 keVee: 

1. A uniform background with free normalization. 

2. Gaussian peaks centered at 1.1 keVee and 1.29 keVee, originating from known radioac- 
tive decays, each of which is modeled as a Gaussian with free normalization and width 



given by the detector energy resolution 5E ~ ^(0.0697)2 + 0.000977 £ fl keV from 

inj. 



3. A contribution from dark matter scattering calculated according to (35), using the 



time-average earth velocity and the quenching factor (37). We do not include any pos- 
sible channeling in Ge. We neglect the finite energy resolution as it is much narrower 
than the intrinsic energy distribution of dark matter recoils. 
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All expectations are multiplied by the detection efficiency plotted in [7T] as a function of 
energy. As the number of events in each bin is small, we use Poisson log-likelihood rather 
than x 2 to quantify goodness of fit. It is important to remember that the CoGeNT low- 
energy excess is entirely consistent with an exponentially falling background. However, we 
do not find it useful to include such a component when fitting the spectrum because it is not 
distinguishable from a dark matter signal with the present experimental statistics. At higher 
statistics, a spectral feature (which is expected in much of the down-scattering parameter 
space) could help to establish the dark matter hypothesis if it is correct. 

DAMA/Nal and DAM A/LIBRA The DAMA/Nal [23] and DAMA/LIBRA [Ml 125] ex- 
periments have reported an annual modulation in the rate of single-hit events with ionization 
energy in the interval 2-5 keVee, with a consistent phase over 13 years. Here we fit this 
signal assuming it originates from dark matter scattering off the Sodium and Iodine nuclei 
in their detectors. 

The detector response to nuclear recoils is complicated by the "channeling effect" [72J, 
whereby a fraction of recoiling nuclei deposit almost all of of their energy electronically. 
Following the parametrization of [73J, we take E ee = En in a fraction / c han of scattering 



events (the "channeled events"), and (37) in the remainder, with 



e -E R /l8 e -E R /40 

/chan,Na^ : - n 7 - = , fchan,!^ , , n Rt ; p > QF(Na) = 0.3, and QF(I) = 0.09. (38) 

Thus, some fraction of events at 2— 5 keVee arise from 2—5 keV nuclear recoils. In the models 
we consider, channeled I recoils dominate the expected signal, with a smaller contribution 
from un-channeled Na recoils. 



It is important to note that the functional forms in (38) are parametrizations of the Monte 
Carlo results in [49] , which differ significantly from analytical models [711 E5] and have not 
been tested experimentally in the relevant energy range. Moreover, partial channeling is 
not included in these parametrizations. All of our quantitative predictions for the DAMA 
experiment depend very sensitively on the behavior of channeling in the 2-5 keV range. 

We calculate the modulation spectrum by subtracting the expected energy spectrum 
when the Earth's velocity in the galactic frame is at its annual minimum, from the expected 
spectrum when this velocity is at its annual maximum. We smear this spectrum with a 
gaussian whose width corresponds to the DAMA energy resolution, 5E = 0.0091 + 0.448/ yE 
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XENON10 The XENON10 experiment used a liquid Xenon target with an exposure of 
316 kg-days, and saw no events consistent with light dark matter. We use the 2-20 keV 
nuclear recoil energy bin and the reported software cut acceptances and nuclear recoil band 
acceptance in [S] to set a limit. 

A large source of uncertainty comes from the scintillation efficiency, L e s, which corre- 
sponds to the amount of nuclear recoil energy that gets recorded as scintillation light. This 
affects the detector's threshold energy, and a small change in L c q has a large impact on the 
resulting constraints. Unless otherwise stated, we use the recent measurement of L c r in [45J , 
which yields weaker constraints than the value L e g = 0.19 quoted in [H]. Our fit to this 
L e ff curve implies a true energy threshold of 4.2 keV r , and = 32.2 MeV. 

CDMS-Si The CDMS experiment contains both Germanium and Silicon detectors. Far 
less data has been collected for the Silicon detectors, but we have checked that this data 
nonetheless sets a tighter constraint than the recent Germanium analysis |77j . For CDMS-Si, 
a total exposure of 53.5 kg-days has been obtained with a threshold of 7 keV nuclear recoil 
[17] (we use this updated result rather than the results in [7S]), corresponding to = 19.2 
MeV. We use the efficiency provided in [47] . While the resulting constraints generate tension 
with the CoGeNT and DAMA preferred regions, they are extremely sensitive to the precise 
threshold energy, as we will demonstrate. 

D. Results 

Figure [T] shows the regions consistent with the CoGeNT and DAMA data for different 
choices of the mass splitting 5 from to 35 keV. We take f p = l and = 0, as appropriate 
for dark matter interacting with nuclei through a kinetically mixed mediator. The elastic 
scattering regions (green region and contours) are consistent with the results for f p = f n = l 
in [261 E7J, after rescaling our results by a factor of (Z(Ge)/A(Ge)) 2 ~ 0.19. The CoGeNT 
regions correspond to regions where — 2Aln£ + 21n£ max = 6.25 (and 14.16 for the right 
plot), where C is the Poisson likelihood and L is its maximum value for any \8\ (here near 
5 = 0). Likewise the DAMA contours correspond to x 2 ~ xL = 6-25 (14.16) relative to 
the overall minimum x 2 , occurring near \S\ = 20 keV. These boundaries approximate formal 
90% (1.6a) and 99.7% (3<r) likelihood regions, but they neglect systematic uncertainties in 
the dark matter halo and detector response that are known to be important. 
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We note that the down-scattering parameter space allows simultaneous fits to the Co- 
GeNT and DAMA data. In our parametrization of channeling, the preferred region is 



m dm -4-10 GeV, a Q x / ~ (1 - 2) X 10~ 40 cm 2 , \8\ -15-35 keV. (39) 

but would shift according to any change in the treatment of channeling. 

Figure [2] shows expected spectra at CoGeNT and DAMA for a benchmark point, \5\ = 25 
keV, mdm = 5.0 GeV, and a ■ f ~ 1.4 x 10~ 40 cm 2 , a point in this preferred region (indicated 
by the dot in Figure [T]). The CoGeNT recoil spectrum agrees with (33) after accounting 



for the quenching factor. The modulation at DAMA comes from the high-energy tail of 
the nuclear recoil energy distribution, and is dominated by channeled iodine recoils (with a 
smaller contribution from unchanneled sodium recoils). The large component below the 2-5 
keVee DAMA region of interest (and below the 2 keV threshold of the DAMA modulation 
analysis) comes from unchanneled iodine scatters. For these parameters, we expect about 4 
events in CDMS-Si and none in XENON10 under our baseline assumptions. Under different 
assumptions about the DAMA detector response, the preferred regions move (and can reduce 
the tension with CDMS-Si and XENON10). 

Constraints CDMS-Si, XENON10, and the DAMA unmodulated single-hit spectrum 
all have the potential to constrain these models. However, significant systematic uncertain- 
ties apply to each of these limits. Curves corresponding to 2.3 and 5 expected events in 
XENON10 and CDMS-Si data, with the baseline assumptions detailed above, are included 
in the right panel of Figure [Tj These limits, especially the limit from CDMS-Si, exhibit ten- 
sion with the preferred overlapping region of DAMA and CoGeNT. However, their severity is 
rather sensitive to various assumptions about each detector's response near threshold, while 
the precise regions preferred by CoGeNT and DAMA are vulnerable to similar uncertainties 
as well as energy-scale uncertainties (and particularly, in the case of DAMA, to the detailed 
features of channeling). We illustrate the sensitivity of the XENON10 and CDMS-Si limits 
to the energy scales of these experiments in Figure [3j where we have varied the true energy 
scale of their detector thresholds by 20%. 

The potential importance of the DAMA unmodulated single-hit rate [21] (which varies 
slowly between 1 and 1.5 counts/kg/day for energies above 1.5 keVee and rises rapidly at 
lower energies) in constraining models of light dark matter was emphasized in [13]. With 
the standard channeling assumptions, our benchmark points overpredict the total rate near 
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FIG. 2: A fit (solid red line) to the CoGeNT data (left) and DAMA data (right). In the CoGeNT 
fit, we assume an energy-dependent detection efficiency [26J and include four contributions: dark 
matter (red dashed), two gaussians for the cosmogenic peaks near 1.1 keV and 1.29 keV (green 
and blue dashed), and a constant background (gray dashed). The dark matter mass (5.0 GeV), 
cross-section (1.4 x 10~ 40 cm 2 ), and splitting 5 = 25 keV correspond to a benchmark point that 
fits both data sets very well (indicated by the dot in Figure [TJ . The expected number of events in 
CDMS-Si and XENON10 is about 4.3 and 0, respectively. Assuming a 20% larger energy threshold 
in CDMS-Si would lead to 1.1 events expected. 

2 keVee by about a factor of 2. This is a potentially severe source of tension, but it is 
driven by two poorly constrained features of channeling — the behavior of channeling at 
low (< 3 keVee) energies, and a possible energy dependence of the quenching factor in this 
range. Blocking and de-channeling effects can dramatically reduce the channeling fraction 
at low energies, rendering the predicted single-hit rate consistent |74j . A careful study of 
channeling in the few-keV region is needed to make sharp predictions of these rates and 
spectral profiles in order to settle this issue. The single-hit tension is also driven by fitting 
the dip in the lowest DAMA bin. Allowing for a mild uncorrected decrease in efficiency 
in this bin (so that the rate in this bin is comparable to the 2nd bin) would render lower 
splittings consistent with the data, and make the unmodulated spectrum consistent. We 
also note that the lower-mass parameter regions discussed in [29], which must assume Ge 
channeling to explain the CoGeNT excess, do not have this problem. 

Finally, although we have checked that the preferred regions and constraints depend only 
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FIG. 3: The effect of energy thresholds on XENON10 and CDMS-Si Limits. Left: The contours 



and thin solid lines are as in the left panel of Figure [TJ The dashed black and gray lines indicate the 
change in XENON10 and CDMS-Si sensitivity when their threshold energy is increased by 20%. 
In the case of CDMS-Si, we use the same rescaling of recoil energies in determining the detection 
efficiency. Right: The spectra in CDMS-Si for the benchmark point with 8 = 25keV (top panels 
of Figure [2]), assuming the nominal efficiency (blue solid curve) or assuming a 20% offset of the 
energy threshold and energy-dependent efficiency (dashed green curve). The total scatter rate with 
no efficiency correction applied (thin dashed gray curve) is also shown for reference. 

mildly on the mean halo velocity, and are quite insensitive to variations in t> esc , other changes 
to the halo velocity profile may have a more significant impact. 

VI. ^-FACTORY PREDICTIONS 

For all of the scenarios outlined above, the mass of the vector A' and higgs states are 
underneath the energy threshold for direct production at the S-factories BaBar and Belle 
[3U |33]. The main production processes are e + e~ — > A /r y and e + e~ — > A'h^. For ~ 
1 GeV, the dominant decay modes of A' are e + e~, or 7r + 7r~, all of which can have 

branching fractions > 10%. Other decays with less than < 10% branching fractions include, 
tt°uj, K + K~, 7r + 7r~27r°, 27r + 27r~, 7r°7r + 7r~, and K°K°. The hd will mainly decay to an A 1 
and two soft leptons or pions, or sometimes directly to two muons or pions [34, 35J. Other 
hadronic final states are also possible for rrih d ^ 1 GeV. From [3U|35], the production cross- 
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FIG. 4: Constraints derived from the e + e~ — > A'j process in the /i + /u~7 search in BaBar's T(3S) 
data using [30]. The gray region is excluded at more than 90% confidence, while the dotted line 
shows the combined reach of BaBar and Belle. The blue region shows the expected range of 
parameters in models of light sub-dominant higgsino dark matter. While not shown, we expect 
that multi-lepton searches sensitive to the reaction e + e~ — > A'hj in the four or more lepton/pion 
final state can reach sensitivities of order (ecw) 2 ~ 0(1CT 8 ) [79J. 



sections are, 



cm 

2^2 \ / 77i2 i ^ \2 



valid over the range 1 GeV < m^/ + m hd < E cm . See [SJ ES] for a detailed discussion of 
these reactions. 

Assuming that the branching fraction for the decay A' — >■ is 1/(1 + < 30% 

for rriA' ^ 1 GeV, current limits on e + e~ — > A'j — > require (ecw) 2 ^ 2 x 1CT 5 over 

the range 1 GeV < vn^ < 6 GeV [40J , and are somewhat stronger at higher mass. Figure 
^illustrates the limit on (ecw) 2 that we derived from . Limits on e + e~ — > 41 JTD] and 
e + e~ — > invisible + 7 (EDI E3] can also be constraining, but are slightly weaker than the 
di-lepton channel for the range of parameters we are considering. Likewise, rare decay 
constraints are also not as strong as the di-lepton search [38, 82J. 
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In the models considered in this paper, (ecw) 2 is predicted as a function of an and m^. 
An improved search for e + e~ — > A'j in the /i + /z~7 final state using the full T(4S*) data set 
at BaBar should be sufficient to probe almost the entire parameter space of this class of 
theories. Even better is a combined analysis of BaBar and Belle data in this channel. If 
kinematically accessible, the e + e~ — > A'hj process can offer a particularly striking signature 
in the six lepton (and pion) final state. Existing searches in four lepton channels [79J indicate 
that upcoming searches in higher lepton (and pion) multiplicity channels should thoroughly 
explore this region. 

For the scenario where the relic abundance of higgsino GeV-Matter arises from thermal 
freeze-out, normalizing «b a m 2 A ,/m~ h to the CoGeNT signal implies a particularly sharp 



prediction for (ecw) 2 , given by (21). The predicted 5-factory cross-sections in this case are 
remarkably close to current limits jlQl I79H8T] and are readily tested by extending existing 
5-factory searches. 

While we believe that 5-factories are uniquely suited to explore the light dark matter 
scenarios outlined in this paper, LEP and hadronic collider searches for lepton jets can also 
play an important role. If Standard Model super-partners are light enough to be directly 
produced at colliders, then supersymmetric dark sectors can give rise to spectacular lepton 
jet signals at the Tevatron or LHC [831487] . However, note that the super-partner scale can 
be as high as ~ m^j^fe ~ 10 TeV before irreducible super-symmetry breaking feeds into the 
dark sectors we consider. Additional searches in J/\l/-factories [H1I12], 0-factories [3T| 157], 
fixed-target experiments [BHHSS], and satellite searches [HH [SU EZ] can provide evidence for 
(or constrain) GeV-scale gauge sectors coupled through kinetic mixing. 



VII. CONCLUSION 



If supersymmetry stabilizes the weak scale, then kinetic mixing with a new gauge force 
can naturally generate a GeV-scale sector. This remarkable possibility could open the door 
to exploring supersymmetry at low energies. In simple supersymmetric dark sectors, stable 
GeV-Matter can arise in the dark higgs sector, and these may comprise a portion of the 
dark matter with new and interesting phenomenology. This possibility is well-motivated by 
cosmic ray data suggestive of GeV-scale forces. 

In this paper, we have proposed simple models of GeV-Matter within a supersymmetric 
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dark sector that couples to the Standard Model through kinetic mixing. The low-energy 
effective Lagrangian we have studied has already appeared in the literature to realize inelastic 
up-scattering of (9(100 GeV) mass dark matter as an explanation of the DAMA signal. Our 
interpretation of the data is different: 

• The CoGeNT and DAMA direct detection anomalies are simultaneously explained 
by inelastic down- scattering of 0( GeV) mass dark-sector states (GeV-Matter), for 
which the higgsinos of a supersymmetric dark sector are an excellent candidate. These 
states comprise a very sub-dominant fraction of the dark matter. Tension with existing 
experiments is reduced, and the CoGeNT and DAMA rates readily agree — in contrast 
to elastic scattering explanations. 

• Electroweak symmetry breaking triggers breaking of the dark sector U(l)u, giving the 
vector multiplet a mass tha 1 ~ ^J^Qd^w ~ GeV. The direct detection cross-section 



of GeV-Matter is predicted (equation (10)) in terms of Standard Model parameters 
with no dependence on e or a^. 

The correct CoGeNT and DAMA event rates are obtained for GeV-Matter comprising 
0.2-1% of the halo mass density. This abundance can be explained by thermal freeze- 
out, and the couplings required for this scenario are readily tested with existing B- 
factory data. 



• Another possibility is that a metastable WIMP decays to GeV-Matter after thermal 
GeV-Matter annihilation has frozen out. In this case, Pllght ~ m '' ght robustly pro- 

PWIMP fWIMP J 1 

duces the desired relic abundance for metastable WIMP masses of 2-3 TeV. We have 
exhibited a model where the metastable WIMP is the bosonic superpartner of the 
TeV-mass dark matter whose decays can explain the cosmic ray data. 

Thus, in our unified and predictive framework, the GeV-scale dark sector plays a crucial 
role in explaining cosmic ray excesses and contains stable GeV-Matter responsible for direct 
detection signals. The most important test of our interpretation of DAMA and CoGeNT 
will come from further searches with direct detection experiments. In particular, further 
studies using CDMS Silicon and Germanium data, especially with lower threshold energy, 
will be vital. The power of these experiments depends sensitively on their threshold ener- 
gies and signal efficiencies. Future studies should include these uncertainties when quoting 
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constraints. Low-energy flavor factories play an equally important role in testing our pro- 
posal. A resonance search in e + e~ — > 7/i + /i~ in the full BaBar or Belle datasets and higher 
lepton-multiplicity searches should be powerful enough to discover or exclude the light vec- 
tor A' over most of the parameter space that can explain CoGeNT and DAMA. Additional 
searches in J/\l/-factories [UJ, ^-factories [36j[37], and fixed-target experiments [88HM] can 
provide evidence for (or constrain) a GeV-scale gauge sector coupled through kinetic mix- 
ing. Finally, if Standard Model super-partners are light enough to be directly produced at 
colliders, then supersymmetric dark sectors can give rise to spectacular lepton jet signals at 
the Tevatron or LHC [SMS EH]- 
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Appendix A: A UV completion for GeV-Matter from Decay 



For completeness, we present a simple UV completion of the model of Section IV, where 
the GeV-Matter states are populated from the decay of a TeV scale state. Consider the 
following matter content: 
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where ip is a real number chosen in such a way to guarantee the form of the superpotential 
and Kahler potential corrections. At some very high scale, the complete superpotential is 
given by 

W = W MS SM + Wdark + W mix + W sp]it + W K + lU break , (Al) 
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where M^mssm is the supersymmetric Standard Model, 

W dark = ~Wl + \Sh u h d + M$A (A2) 

contains the GeV-scale dark sector fields (including the higgsinos, which will form the GeV- 
Matter) and the TeV-scale dominant dark matter component from the $ U) d, and 

W mix = - e -W$W aD (A3) 

contains the gauge-kinetic mixing term that connects the supersymmetric Standard Model 
to the dark sector. In addition, 

HVt = Z x h u h d + Z 2 $ u h d - \m x Z\ - X -M 2 Z\ (A4) 

will, after integrating out the ~TeV-mass fields Z\ and Z2, generate the higher dimensional 



terms in Eq. (26) that generate the mass splittings among the heavy and light dark matter 



states. The superpotential 

W K = AX$ U + AYh u + M*AA (A5) 
will, after integrating out the fields A and A at some high scale M*, generate the Kahler 



potential term Eq. (27), but no additional superpotential terms. Note that the symmetries 



guarantee that the Kahler potential has no lower-dimensional terms, or other terms of the 



same dimension as Eq. (27). Finally, we will not specify the precise form of WVeak- This 
superpotential contains fields that break supersymmetry to give non-zero F-terms Fx and 
Fy- It also has to contain the supersymmetry breaking sector that generates the soft masses 
in the supersymmetric Standard Model via some mediation mechanism. 

Appendix B: Sources of WIMP and GeV-Matter Mass Splittings 

We briefly remark on three models that give rise to inelastic splittings for the $ fermions 
and dark higgsinos. The simplest model has charges ±1/2, so that additional interac- 
tions 

W = \+h d $l + \-h u $ 2 d (Bl) 

are possible. The former directly generates a somewhat large Majorana mass for the heavy 
dark matter, while integrating out $ loops or other heavy states generates an operator of 



the form (16) and therefore a splitting of the dark higgsinos. 
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Smaller hierarchies of splittings are achieved in models where $> U) d have charge ±1. We 
first consider a model with two singlets X and Y, and two Z 2 global symmetries: one under 
which X, S, and h u are odd and one under which Y, S, and hd are odd. These symmetries 
allow, in addition to ^ and (11), relevant and marginal interactions 



W XY = HxX 2 + fi Y Y 2 + \^ d h u X + \ + $ u h d Y + X S SXY. (B2) 

Integrating out X and Y gives rise to interactions 

— &hl and — hlhl, (B3) 

which split the heavy dark matter and dark higgsinos respectively. The $ splitting is gen- 
erated at tree-level and the dark higgsino splitting only at loop-level, so the dark higgsino 
splitting is naturally smaller. For Yukawa couplings ~ 1/10, (hd) ~ 10 GeV, and M ~ 1 
TeV the typical splittings are 1 MeV and 10 keV respectively for the heavy and light states. 

Alternately, our $ and U(1)d higgs-sector fields can couple to a doubly-charged field T±2 
through interactions 

W T = fi T T +2 T. 2 + X d h 2 d T +2 + \ u h 2 u T_ 2 + A_$^T +2 + A + $2T_ 2 . (B4) 

After U(1)d symmetry-breaking, the F-flat condition for T +2 is satisfied if T_ 2 acquires a 
small vev (T_ 2 ) = —XdV^/^ij- ~ MeV for fi? ~ few TeV. This vev induces small Majorana 
masses for both $^ and h d ] a hierarchy of Yukawa couplings X u A + leads to a hierarchy 
of splittings. 



Appendix C: Form Factor and Halo Parametrizations for Direct Detection 



We elaborate here on the form factor and halo parametrizations appearing in (35), which 



are used throughout Section |Vj F(Er) is a nuclear form factor, which we model as a Helm 
form factor as in |100j . 



F{qr r , 



3ji(gr n ) ^_ (9s) 2 /2 



qr ri 



(CI) 



with q = y/2m N E R , r 2 = c 2 + |7r 2 a 2 
fm. 



5s 2 , c = 1.23A 1 / 3 - 0.60 fm, s = 0.9 fm, and a = 0.52 



In the final factor of (35) 



*w f(y-y e ) 

a v . 



(C2) 
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v is the dark matter velocity in the frame of the Earth, v e is the velocity of the Earth in 
galactic coordinates, and f(v — v e ) is the local dark matter velocity distribution in the Earth 
rest frame, which we model by 



/( 



v — V P oc 



-(v-v e ) 2 /vl 



(C3) 



for v < v esc and otherwise. We normalize f(v — v e ) to unity (for example, if t> esc = oo, 
then the proportionality constant is 1 / (7TVq) 3 ' 2 ) . We use vq = 254 km/s |101j for both the 
average velocity of local stars and the mean velocity of halo particles. The local galactic 
escape velocity v esc is measured to lie in the range 498 km/s < v esc < 608 km/s at 90% 
confidence |102j . and we use the median likelihood value of v esc ~ 544 km/s. The Earth 
velocity in the galactic frame is determined as a function of time by v e = v Q +t>©, where v & 
is the Sun's velocity relative to the Galactic rest frame |1011 1103] , 



( 10.00^ 



\ 



5.23 
7.17 



km/s + 



V / 



(C4) 



and Vq is the Earth's velocity relative to the Sun, 

/ 

% = (u E )(l - esin(A - A )) 



km/s. 



(C5) 



cos ((3 X ) sin (A - X x ) 
cos((3y) sin(A — X y ) 
\ cos(/3 2 ) sin(A - A 2 ) I 

Here, the Earth's orbit has a mean velocity (ue) = 29.79 km/s and ellipticity e = 0.016722, 
A is the angular position of the Earth's orbit, and the quantities f5 x ,y,z and \o, x ,y,z are given 
in [TOO] . 

The velocity integral in (35) cuts off at a minimum velocity v m i n determined by the 
scattering kinematics: 



m N E R 



+ 5 



(C6) 



y/2m N E R Hn 

where fijy is the dark matter-nucleus reduced mass, and 5 is the splitting between the 
incoming and outgoing dark matter particles: 5 = for elastic scattering, 5 > if the 
incoming dark matter particle scatters into a heavier excited state ("up-scattering") [104J , 
while 5 < if the incoming dark matter particle is already in an excited state and scatters 
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into a lighter ground state ("down-scattering"). For the masses and splittings we consider, 
at most a tiny fraction of the halo can up-scatter and the down-scattering signal dominates. 
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